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Abstract. Equilibrium thermodynamic laws are typically applied to horizons in 
general relativity without stating the conditions that bring them into equilibrium. 
We fill this gap by applying a new thermodynamic interpretation to a generalized 
Raychaudhuri equation for a worldsheet orthogonal to a closed spacelike 2-surface, 
the “screen”, which encompasses a system of arbitrary size in nonequilibrium with its 
surroundings in general. In the case of spherical symmetry this enables us to identify 
quasilocal thermodynamic potentials directly related to standard quasilocal energy 
definitions. Quasilocal thermodynamic equilibrium is defined by minimizing the mean 
extrinsic curvature of the screen. Moreover, without any direct reference to surface 
gravity, we find that the system comes into quasilocal thermodynamic equilibrium when 
the screen is located at a generalized apparent horizon. Examples of the Schwarzschild, 
Friedmann-Lemaitre and Lemaitre-Tolman geometries are investigated and compared. 
Conditions for the quasilocal thermodynamic and hydrodynamic equilibrium states to 
coincide are also discussed, and a quasilocal virial relation is suggested as a potential 
application of this approach. 


PACS numbers: 04.20.Cv, 04.70.-s, 04.70.Dy, 98.80.-k 


Class. Quantum Grav. 32 (2015) 165011 


Quasilocal energy and thermodynamic equilibrium conditions 2 

1. Introduction 

A thermodynamic description of general relativity has been a long-sought goal CD 13 
which intensified with the advent of black hole mechanics M. Most studies in 
the literature focus on equilibrium thermodynamic^ of horizons , without stating the 
conditions that bring them into equilibrium. In fact, in gravitational physics there 
are no well-defined conditions for defining equilibrium in terms of the behaviour of a 
system itself. In this paper we will take steps to remedy this by defining a quasilocal 
thermodynamic equilibrium condition using a purely geometric approach. We will focus 
on the extrinsic geometry of the closed spacelike 2-surface that appears both in various 
quasilocal energy definitions [T0HI2] and the generalized Raychaudhuri equation of 
Capovilla and Guven [13]. This connection provides a natural way of defining quasilocal 
thermodynamic potentials, and hence equilibrium conditions. 

In general, if one wants to investigate the energy exchange mechanisms of a 
gravitational system from the thermodynamic viewpoint, the system should have a 
finite spatial size. Energy definitions which refer to the spatial asymptotic behaviour 
are not good candidates for general thermodynamic equations. Thus quasilocal energy 
definitions, which refer to a Hamiltonian on the 2-dimensional spacelike boundary 
mm, are very important for general relativistic thermodynamics. Here we will 
link such definitions to a generalized notion of the work done in the deviation of 
worldsheet congruences, to define quasilocal thermodynamic potentials. We also provide 
a quasilocal first law by considering a worldsheet total variation, in which a quasilocal 
temperature can be understood as a worldsheet-constant. 

We will consider general spherically symmetric spacetimes, rather than just black 
holes. Naturally, the thermodynamics of more general spacetimes is also of interest. 
However, quasilocal definitions require closed surfaces on which to integrate to define 
the system in question. Thus we believe that the new approach taken here could 
also be applied in future to systems whose bounding surfaces are defined by different 
symmetries. 

Although early investigations dealt with equilibrium thermodynamics of black hole 
event horizons, in the last few decades more general trapping, apparent and dynamical 
horizons of generic spacetimes have been introduced HM9] . Hayward’s construction 
of equilibrium thermodynamics on trapping horizons |15] highlights the significance of 
generalized apparent horizons. While the original definition of apparent horizons applies 
to black holes which require asymptotically flat spatial hypersurfaces [SO], Hayward’s 
generalized apparent horizon which was first constructed for black holes, has also been 
applied in more general cases P2U2HE2!. These include cosmological applications, where 
the generalized apparent horizon is not necessarily spacelike but can be timelike or null 
depending on the equation of state of the cosmic fluid [fl9] . 

In this paper we will consider a spherically symmetric gravitational system of 
arbitrary size which is not in equilibrium with its surroundings. As one of our results we 


t See [6H9j for some exceptions. 
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will show that when a particular equilibrium condition is applied to such a system then 
the 2-surface enclosing the system is located at the generalized apparent horizon of [15]. 
This result makes no direct reference to the surface gravity, which is conventionally used 
to define the temperature of the horizon. 

The paper is constructed as follows. In section [2] we review the quasilocal energy 
definitions in literature, highlighting those features important for our new construction, 
including in particular the determination of a well-defined internal energy. Section [3] 
starts with a short summary of the geometric formalism developed by Capovilla and 
Guven pL3] and continues with interpretation of their generalized Raychaudhuri equation 
as a thermodynamic relation. Following this, a quasilocal thermodynamic equilibrium 
condition and the corresponding thermodynamic potentials are introduced. In section 
[4] these results are applied to the Schwarzschild, Friedmann-Lemaitre-Robertson-Walker 
and Lemaitre-Tolman spacetimes. In section [5] we highlight the difference between local 
thermodynamics of matter fields on curved background and quasilocal gravitational 
thermodynamics, as a precursor to suggesting a potential application of our approach 
to a quasilocal virial relation in section [6] We will use natural units in which c, G,h,ks 
are taken to be 1 throughout the paper. The metric signature will be (—, +, +, +). 

2. Quasilocal energy definitions in literature 

Here, we will discuss those quasilocal energy definitions which will be most relevant 
for our investigation of equilibrium thermodynamic conditions. A detailed review of 
quasilocal energy definitions can be found in [ 23] . 

2.1. Misner-Sharp-Hernandez (MSH) Energy 

In horizon thermodynamics, there is a broad consensus [T9II24] [27] on the choice of the 
internal energy of a generic spherically symmetric spacetime. It is usually taken as the 
Misner-Sharp-Hernandez energy [28[[29] which we now define. Consider a spherically 
symmetric spacetime metric with coordinates {x a ,9,(p} where {x a } = {f, r}, 

ds 2 = 7 a gdx a dx^ + R 2 (x) (d6 2 + sin 2 Odcj) 2 ) , (1) 

R being the areal radius. In order to study time evolution, one can pick a preferred 
timelike vector, called the Kodarna vector [30], which can be used to define surface 
gravity for dynamic spherically symmetric spacetimes cat Surface gravity, up to a 
constant, is in general related to the temperature defined on the horizon. The Kodarna 
vector is unique and it is parallel to the timelike Killing vector in static spacetimes. Its 
components are given by, 

K a (x) = e aP dyR, K° = 0, K* = 0, (2) 

where e a g is the Levi-Civita tensor in 2-dimensions. Now consider a 3-dimensional 
spacelike hypersurface, £, with induced coordinates, {y a }, induced metric, h aK , and 


Quasilocal energy and thermodynamic equilibrium conditions 


4 


unit normal, n", aligned with the Kodama vector. The Kodarna vector is associated 
with conserved charges including an energy 

E M sh = J cPyy/hT^KW, (3) 

where is the stress-energy tensor of matter in the 4-dimensional spacetime. This 
defines the Misner-Sharp-Hernandez energy which can also be written 

Emsh = | (1 - 7 af> d a RdgR) . (4) 

2.2. Brown-York (BY) Energy 

Brown and York followed a Hamilton-Jacobi approach to define a quasilocal energy 
on a 2-dimensional spacelike boundary, §< = <9£ t , of a compact spacelike 3-dimensional 
hypersurface, Y t . These hypersurfaces are taken to foliate a compact spacetime domain, 
D, with topology £ x The fundamental object of the Brown-York definition is 

the 3-dimensional timelike boundary of the domain D which will be denoted by B. 
The boundary B is foliated by § t and it admits a fixed metric. According to this 
approach, the stress-energy tensor defined on B captures the coupled effects of matter 
and gravitation. When one projects this stress-energy tensor tangentially and normally 
to the spacclike 2-boundary of B, one obtains the quasilocal energy, momentum and 
spatial stress density. The Brown-York quasilocal energy density (energy per 2-surface 
area) can be written as the following 

e = (k - k 0 ), (5) 

where k is the extrinsic curvature of § when it is embedded in £ and k$ is the 
corresponding extrinsic curvature evaluated for a suitable reference spacetime. Starting 
with the action principle, one can show that the Brown-York quasilocal energy is the 
value of the 2-dimensional boundary Hamiltonian once the constraints are introduced 
in the bulk. The Brown-York energy is given by 

E by = (£d§>(k - k 0 ). (6) 

Svr 


2.3. Kijowski (K) Energies 

Kijowski [TO] followed a Hamilton-Jacobi approach which is fundamentally different 
to that of Brown and York. He attacked the quasilocal energy problem from 
first principles, by considering a symplectic structure on the phase space which is 
consistent with a nonvanishing boundary HamiltoniaiJ§|. This can be interpreted as 
a modification of the Arnowitt-Deser-Misner symplectic structure which assumes the 
vanishing of the boundary Hamiltonian under certain asymptotic conditions. In fact, 


§ See ]3l] for a similar viewpoint. 
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in collaboration with Tulczyjew, Kijowski used these symplectic structures to define 
boundary Hamiltonians [221133] decades before other quasilocal energy definitions were 
introduced. This work appears to have been neglected by the community until Kijowski 
applied his approach to general relativity explicitly [XU] . 

Kijowski applied Legendre transformations to the boundary Hamiltonian and 
obtained two physically meaningful quasilocal energy definitions. He called them the 
internal energy and free energy due to them being related to each other via Legendre 
transformations with Dirichlet and Neumann type boundary conditions. The two 
quasilocal energy definitions of Kijowski [ID] are respectively, 


rpin _ 

h K - 


I67r 


(e - p) - ki 


( 7 ) 


E k — — - 


8tt 




( 8 ) 


where k and l are the trace of the extrinsic curvature of § with respect to its spacelike 
and timelike unit normals. Thus \J(k 2 — l 2 ) gives twice the mean extrinsic curvature of 


One may calculate ko by isometrically embedding § into Minkowski spacetime 
hyperplanes and dragging the 2-surface along the timelike Killing vector of the 
Minkowski spacetime. The reference energy density, ko, is the extrinsic curvature of 
the worldtube obtained in this manner and it is related to the boundary conditions. 
Note that E l ff reduces to the Misner-Sharp-Hernandez energy for spherically symmetric 
spacetimes and gives the Hawking mass for the Schwarzschild geometry [23] . Therefore 
E l ff will be taken as the internal energy of the quasilocally defined system in this paper. 
However, we do not follow Kijowski in referring to Ek as the free energy. This is 
because it takes the same form as the Liu-Yau quasilocal energy HU which is interpreted 
differently in the literature, as we will now discuss. 


2-4- Liu-Yau (LY) Energy 

In Liu and Yau’s work HU there is no reference to a timclike 3-dimensional boundary. 
They considered the embedding of § directly into a spacetime domain D by taking its 
two normal null vectors and the corresponding mean extrinsic curvature. That provides 
a well-defined quasilocal energy under the direct embedding of a 2-dimensional spacclike 
surface into a 4-dimensional spacetime. This is the motivation for the method we use. 
When Liu and Yau converted from their original notation into the one used here, their 
energy expression becomes: 


Ely = ~ 


1 

87T 


VW - p) - fc) 


( 9 ) 


The reference energy is obtained by embedding § into the 3-dimensional Euclidean space, 
M 3 , and calculating its extrinsic curvature, ko- This isometric embedding is unique up to 
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the isometries of M 3 . Note that their quasilocal energy expression is equal to Kijowski’s 
energy given by (jHJ). The positivity of the Kijowski-Liu-Yau energy, denoted E KLY from 
now on, has been proven by Liu and Yau mm- It is a widely accepted quasilocal 
energy definition [231136] . 

2.5. On radial boost invariance 

The quasilocal energies Eff and Ekly of a system are obtained via the mean extrinsic 
curvature of a 2-dimensional spacelike boundary, which we will call the screen. For a 
spherically symmetric spacetime, for example, the generator of the quasilocal energy is 
not just a single timelike vector. Rather, one needs to consider both the future pointing 
timelike normal and the outward pointing radial spacelike normal in order to calculate 
the mean extrinsic curvature of the screen. 

Following the ideas of Epifiida 

Ekly can be interpreted as a proper mass-energy 
of the system. It is known that both Ekly and Eft are invariant under radial boosts 
of the quasilocal observers who define the 2-surface [T0TH21 [23 j. Such a property is 
necessary to define a system consistently since one needs the ability to keep constant 
the degrees of freedom that define the screen enclosing the system. In the case here, 
{9, 0} are the coordinates on the screen that are kept constant. Then the evolution 
of the system is investigated by perturbing the screen along the remaining degrees of 
freedom, parametrized by the { t , r} coordinates. Thus the screen observers agree on 
the quasilocal energy content of the same system irrespective of them being boosted or 
having instantaneous radial accelerations with respect to any other screen. 

We note that in this picture, the Kodarna vector is an object that lives in the 
temporal-radial plane. The fact that the Kodarna vector is associated with a conserved 
Misner-Sharp-Hernandez energy has previously been described as a miracle [37]. Here 
we emphasize that since E l £ matches Emsh under spherical symmetry, this miracle is 
a natural consequence of any consistent quasilocal Hamiltonian formalism of general 
relativity. Thus one does not need to define a single preferred observer in the energy 
calculations. 

The distinction of the degrees of freedom that are used to define the system and 
the ones used in the investigation of the evolution is crucial for the interpretation of the 
formalism introduced in the next section. 

|| Note that one can find a topologically spherical 2-surface in flat spacetime with negative Brown-York 
energy and positive Kijowski and Kijowski-Liu-Yau energies. According to O Murchadha, Szabados 
and Tod [35] this might take place if § does not he in a spacelike hyperplane of flat spacetime. However, 
in the present paper we only consider spherical symmetry, in which case the Kijowski and Kijowski- 
Liu-Yau energies are well defined. 

1 Epp defined a similar energy to Ekly with a different nonunique choice of reference energy. 
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Capovilla and Guven ra generalize the Raychaudhuri equation which gives the 
focusing of an arbitrary dimensional timelike worldsheet that is embedded in an 
arbitrary dimensional spacetime. We apply their formalism to a 2-dimensional timelike 
worldsheet, T, embedded in a 4-dimensional spacetime, using their original notation. 

Consider an embedding of an oriented worldsheet with an induced metric, rj ab , 
written in terms of orthonormal basis tangent vectors, {E a }, 

9( E ai E b) =Vab, (10) 

where g lw is the 4-dimensional spacetime metric. Now consider the two unit normal 
vectors, {n*}, of the worldsheet which are defined up to a local rotation by, 

9(n i ,n J ) = S ij ( 11 ) 

g(n‘,EJ = 0, (12) 

where {a,b} = {0,1} and {i,j} = {2,3} are the diad indices and the Greek indices 
refer to 4-dimensional spacetime coordinates. Also note that to raise (or lower) the 
indices of tangential and normal diad indices one should use g ab (or g ab ) and <W (or S t] ) 
respectively. 

At this point it is important to distinguish the different covariant derivative 
operators [13]. Let the torsionless covariant derivatives defined by the spacetime 
coordinate metric be D and their projection onto the worldsheet be denoted by 
D a = E ,l a . On the worldsheet T, we introduce two covariant derivatives. V a is 
defined with respect to the intrinsic metric and V a is defined on tensors under rotations 
of the normal frame on S. To study the deformations of T and §, the following extrinsic 
variables are introduced [13]. The extrinsic curvature, Ricci rotation coefficients and 
extrinsic twist of T are respectively defined by, 

Kj = -9,AD a E'‘ b )n'' l = Kj, (13) 

labc = V ( D a E “l,) E \ = -~tacb, C 4 ) 

<= = V (£>x‘) (15) 

while the extrinsic curvature, Ricci rotation coefficients and extrinsic twist of § are 
respectively defined by, 


Jj 1 = V (D'E\) n->, 

(16) 

lijk = v ( D i n “i) n \ = -7ay. 

(17) 

sj = v ( D ' E “«) = -sj- 

(18) 


By using those extrinsic variables one can investigate how the orthonormal basis {E a ,n' 1 } 
varies when perturbed on T according to, 

E a E b = lab E c — E ab n ii 


(19) 
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Dy = KjE * + ra.% (20) 

or perturbed on § according to, 

D.E, = S aU E b + J aij n>, (21) 

D i n j = ~J*ij E ° + 7i j k n k . (22) 


Then the generalized Raychandhnri equation which quantifies how much a worldsheet 
focuses is found to be [13J 

V t + = - V‘Kj - J b \jy - KjKy + g (R(E b ,n‘)E a y) , (23) 

where R a g llu is the Riemann tensor of the 4-dimensional spacetime Il3j . and 

g(R(E a , n,)E b , n,) = R^E^n^E^,. (24) 

Note that w bi k transforms as a connection under the rotation of 8 and 

^ b^aij = V \J a ij ~ W bj k Jaki~ W bi k Jakj- (25) 

D, J .. — 7, C J 

0 axj 1 ba czj 


Likewise, S ab l transforms as a connection under the rotation of T such that 


V.AV = 




_ c is CJ — 9 K CJ 

^aci^b ^bci^a 


(26) 


D K , J -7 \K , k 

x ab '1 k ab 


If one contracts (l23li with the orthogonal basis metrics 7] ab and 5^ one finds 

(++) rf% = - (V‘KJ) 

+ g(R(E b ,n‘)E a ,ngyS ij , (27) 

which is the equation to which we will assign a quasilocal thermodynamic interpretation 
to arrive at the key results of this paper. 

In the 3+1 formalism, the Raychaudhuri equation [38] for a congruence of timelike 
worldlines tells us how much the congruence expands or contracts. For a timelike 
congruence, that is not necessarily geodesic, with vanishing shear and vorticity, the 
Raychaudhuri equation reads [}39j 

- ^ = +^© 2 + R^v v - D, KZW), (28) 

where 0 is the expansion scalar, A is the affine parameter, R pi/ is the Ricci tensor and 
t+ are the components of the 4-velocity vector field of an observer. In the 2+2 language, 
however, the generalized Raychaudhuri equation (|2T|) tells us how much a congruence 
of worldsheets, rather than the worldlines, focuses. 

Now consider a general spherically symmetric spacetime. For radially moving 
observers, the extrinsic curvature and the extrinsic twist of T vanishes as well as the 
extrinsic twist of S. Then (1271) reduces to 

- (++) = +v+<f+ - g(R(E b ,n<)E a yX b S tj . 


(29) 
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We will interpret (I29p as a thermodynamic relation for a quasilocally defined system. 
Note that the choice of the orthonormal dyad metrics are r/ a & = (—1,1) and Sij = (1,1) 
for the rest of the paper. 


3.1. Quasilocal thermodynamic equilibrium conditions 

At this point one needs to take care with the definitions of the thermodynamic 
variables under equilibrium and nonequilibrium conditions. In general, equilibrium 
can be seen as a specific state of a system that is ordinarily in nonequilibrium with 
its surroundings. Given specific equilibrium conditions, the equilibrium state acts as an 
attractor to bring the system into a preferably stable state [40]. Moreover, in classical 
thermodynamics, only in the equilibrium case are the existence of thermodynamic 
potentials guaranteed [41 . For such equilibrium states the thermodynamic potentials 
are Lyapunov functioned, and they can be written as linear combinations of each 
other [43]. 

Let us recall the definitions of thermodynamic potentials in classical thermodynam¬ 
ics at the equilibrium state PJ, 

Helmholtz Free Energy : F = U — T S, (30) 

Gibbs Free Energy : Q = F + “W, (31) 

Enthalpy : H = Q + H S = U + , (32) 

where U is the internal energy, and “W represents the work terms which may include 

PV type and other types of work in general. The ‘=’ sign will be used for equations 

that hold only at quasilocal thermodynamic equilibrium from now on. 

Note that Helmholtz free energy is the amount of reversible work done on a system 
in an isothermal process B5f It is one of the thermodynamic variables that can be 
defined both in equilibrium and in nonequilibrium states [46] . Moreover, one way 
of defining the thermodynamic equilibrium is to set the Helmholtz free energy to its 
minimum value [43]. For this reason the Helmholtz free energy provides a physically 
natural means to interpret (129|) thermodynamically. 


3.1.1. Helmholtz free energy density We will take the extrinsic curvature of § as a 
measure of the matter + gravitational Helmholtz free energy density and defineQ 

/“/«=2 (•/„«■/“,). (33) 

since 

/ ee JJ^f a = Vtf—p. (34) 

+ Lyapunov functions are nonnegative functions that have at least one local maxima or minima at a 
point of interest. They are continuous functions with continuous first order derivatives and they vary 
monotonically with the evolution parameter [¥2]. 

* The J a * 7 J ra ?i /4 term appears in the definition of the Hawking [17] and Liu-Yau m mass-energies, 
as the term fip in the notation of these authors. 
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is the object that appears in the quasilocal energy definitions of section [21 Then the 
Helmholtz free energy of the system is obtained one e f is integrated on S, i.e., 


7 — 



(35) 


Since the equilibrium condition is defined in this case by the minimum of the Helmholtz 
free energy, other thermodynamic potentials should be written as linear combinations 
of each other once one sets J = lF m in■ Thus at equilibrium, the Gibbs free energy and 
the internal energy should read 


Q = Frrnn + W, (36) 

< U = R min + ‘TS. (37) 

where ‘W, T and S are to be defined. The Helmholtz free energy density defined by 
(1221) and (CUD is required to be a nonnegative real scalar, and the minimum value it can 
take is zero. This brings us to write the equations above with ‘Emin = 0 as 


<? - ■»' 

U = <TS. 

Then when f — Jf a fa — Vk 2 — l 2 — 0, llTHll becomes 


(38) 

(39) 


(40) 


Recalling the fact that thermodynamic potentials are nonnegative real functions at 
equilibrium, we will force the terms in (1401) to take nonnegative values by taking the 
absolute value of each side before we start making further quasilocal thermodynamic 
interpretations. Thus, 


-V J“. 

v a i 


-g(R(E b ,nQE b , nj ) |. 


(41) 


3.1.2. Work density We will now give a thermodynamic interpretation to the quantity 
on the r.li.s. of (14 1 [) . In the 3+1 formalism, when one considers two observers on 
neighbouring timelike geodesics the deviation of the geodesics determines the relative 
accelerations of the observers. If we consider the spacelikc separation 4-vector, £, 
that connects the neighbouring geodesics, then the components of the relative tidal 
acceleration are given by [39j 

d 2 ^ 

(42) 

where r is the proper time. Thus for a spherically symmetric spacetime we define a 
relative work term that mimics W = F ■ x by 

(d 2 ^\ 

V dr 2 J 


^ = R\ p y v 'r£r 


(43) 









Quasilocal energy and thermodynamic equilibrium conditions 


11 


This relative work term can be interpreted as a measure of energy expended within 
the surface of a body to stretch or contract it under the influence of tidal forces, if we 
assume £ lives on the screen, S. 

Our interpretation is similar to that of Schutz [48] who considered the limits of 
validity of the geodesic deviation equation and calculated the second order contributions. 
He also acknowledged the fact that connecting two geodesics with a separation 
vector is essentially nonlocal. Thus the reason (142|) is valid only for nearly parallel, 
neighbouring geodesics is simply due to observers trying to measure a nonlocal quantity, 
locally. Consequently (l43jl has a more fundamental quasilocal interpretation and the 
| — g(R(E b , nQE b , n^) | term on the r.li.s. of (14TT) might be taken as a measure of work 
density attributed to S. 

To understand this intuitively, consider the analogy of a soap bubble. The work 
done per area to create the surface of a bubble is [44] 

= (44) 

where 7 is the surface tension and dA is a surface area element of the bubble. According 
to classical theory, surface tension arises due to the unbalanced intermolecular forces in 
the bubble. Likewise, according to the analogy formed here, | —g(R(E bl nQE b , n^) | is 
a measure of energy density due to the relative tidal forces that observers experience 
when they move along radial worldlines. This is of course applicable for observers who 
share the same screen §. Therefore at quasilocal thermodynamic equilibrium, we can 
define a general work density, namely a type of surface tension, according to 

u) = ^Ju) a u) a = \J 2 \ -g(R(E b ,ni)E b ,n j )\, ( 45 ) 

so that the amount of corresponding work is given by 

<W = — i TV ■ dS. (46) 

16 vr ,/ § 

3.1.3. Gibbs free energy density In classical thermodynamics, when equilibrium is 
defined by the minimum of the Helmholtz free energy, the Gibbs free energy reads [44] 

g d a„ = jl-dA, ( 47 ) 

for the thermodynamics of surfaces with constant pressure. Following the analogy with 
the surface of a soap bubble, 

^dass = 7 ■ dA 77 dV = j) zv ■ dS, (48) 


Qclass 


/ 


so that 


7 • dA 77 Q = dV 


(49) 
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should hold. This is consistent with (1551) which states that Q = ‘W since the minimum 
Helmholtz free energy is zero according to the equilibrium condition defined here. Tims, 
the l.h.s. of flUD can be taken as a measure of the quasilocal Gibbs free energy density: 



(50) 


from which the Gibbs energy can be obtained by 



(51) 


Note that, in general, the Raychaudhuri equation becomes nonlinear if one wants 
to write it in terms of the energy densities defined here. However, recall that the 
existence of thermodynamic potentials is guaranteed only in the equilibrium case in 
which the potentials can be written linearly in terms of each other. In classical surface 
thermodynamics, the surface tension, pressure gradient across the surface and mean 
curvature of the surface can be related via the Young-Laplace equation {39]. Ideally, in 
order to reach the equilibrium, fluids tend to minimize their surfaces until they have 
zero mean curvature. This is when the the surface tension takes its critical value. In 
the formalism presented here, which is in line with our analogy, this happens when 
f = \Jf a fa = Vk 2 — l 2 = 0, which defines the apparent horizon of a given spacetime. 

Here we use a general apparent horizon S3. defined by the marginal surfaces 
on which at least one of the expansion scalars of the null congruences is zero, i.e., 
$(/)$(n) = 0, where l a ( n a ) is the outward (inward) pointing future-directed normal. Both 
the conditions { 6^ > 0, 6 ( n ) = 0} and {9^ = 0, 9 („) < 0} have previously been used 
to define apparent horizons [3911231122] . Here J a u J a jt in (1331) gives a measure of 9^9^ n y 
Thus when it is equated to zero, one can conclude that at least one of the expansion 
scalars of the incoming or outgoing null congruences converges without knowing which 
one actually does. 

3.1.4■ Internal energy density On introducing the quasilocal energy definitions in 
section [2j we stated that E l £ is a good candidate for the total matter + gravitational 
energy content of a system. It is derived via a Hamilton-Jacobi formalism with Dirichlet 
boundary conditions. According to Kijowski those boundary conditions are associated 
with the true degrees of freedom of the quasilocally defined domain that gives the true 
energy |T0j. When the equilibrium condition is imposed, the internal energy density in 
© can be written as 


u = k 0 . 

Thus the quasilocal internal energy at equilibrium becomes 


(52) 



( 53 ) 
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which should satisfy the equilibrium condition (fT 7 !!) without any PV type termlj. This 
requires that we define a quasilocal entropy and temperature at equilibrium. 

Since the 2-surface § located at the generalized apparent horizon enters naturally, 
we can follow the traditional approach of Bardeen [3] and define the quasilocal 
equilibrium entropy: 

. Area (§) 

5 = A 

where A is a constant which is usually taken to be 4 in gravitational equilibrium 
thermodynamics of horizons. By (l39lh at equilibrium the quasilocal temperature of 
the system is then given by 

U 

T=j (55) 

without any direct reference to surface gravity. 

3.1.5. First law of thermodynamics According to the formalism constructed here, the 
first law should be written as 



5U = 5{PS) = FdS, 


(56) 


since we defined the quasilocal thermodynamic equilibrium via the minimization of 
the Helmholtz free energy which is applicable for isothermal processes. The problem 
with some of the gravitational thermodynamic constructions is that the total variation 
of the internal energy and entropy is performed in specific directions for which the 


first law does not have the dimensions of energv lffl . However, in our framework the 


quasilocal behaviour of the system sets the degrees of freedom with respect to which 
the total variation can be defined. These are the degrees of freedom that live on the 
instantaneously defined timelike surface T. Hence, we will set the total variation to be 
the one on the worldsheet and write 


5U = ^\fv a ‘ll\7 a ‘U. (57) 

Thus the first law should read 

^\/v a UV a U = ‘T (^\/v a SV a S^ , (58) 

with 

6? = = 0, (59) 

U When Kijowski applied the boundary conditions in his original derivation |10j . he chose the induced 
metric components of § to be time independent. This choice holds both for E l £ and for Ekly and it 
is set to prevent the extra inclusion of volume acquisitions in the energy definitions. However, Liu and 
Yau do not set such a condition on Ekly- Moreover, according to the equilibrium conditions set here, 
U contains only the tT S term. Therefore in our case E l £ can be safely used for generic coordinate 
representations of spherically symmetric spacetimes at equilibrium, 
ft For example see fl9] . 
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where {a, b} = {0,1}. Thus the temperature will be a worldsheet-constant rather than 
a constant with respect to some coordinate time. For the examples that are presented 
in the next section, it is easy to check that (j59j) is satisfied. Note that, in general, 

E» a D, -> E\d^ (60) 

can be used for the scalar functions that appear in (1581) and (1591) . 

4. Examples 

4-1. Schwarzschild geometry 

Consider the Schwarzschild metric in standard coordinates 

<is 2 = - (5 - </( 2 + (V + rW, (61) 

where dQ 2 = d9 2 + sin 2 9d(j) 2 . As stated previously, to define a consistent system observers 
should have fixed angular coordinates. As one example consider static radial observers 
with double diad 

^=( 0 ' 0 ' 0 ’^)' (63) 

The choice of static observers is inconsequential for our results, which also apply to 
observers with an arbitrary instantaneous radial boost with respect to this frame. For 
such observers, (l33l) implies 



which can be substituted in the Raychaudhuri equation, (|29lh in the general 
nonequilibrium case to give a notion of nonequilibrium quasilocal energy exchange. 

In order to set the quasilocal thermodynamic equilibrium condition, 4F should be 
minimized. It is easy to see that, this occurs when r = 2 M which coincides with 
the location of the black hole horizon. Now let us calculate the internal energy at 
equilibrium. Given the metric in (16TT) and the isometric embedding of § into Euclidean 
3-space, one finds k 0 = 2/r. Then according to (1391) and (1531) . 

<U = ‘TS = M, (65) 

with T = A/(87rr) = A/(167 tM), and S = Area( §)/A = (167 tM 2 )/A where A is a 
constant. For those who wish to relate this temperature to the Hawking temperature 
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[50], there is a problem of factor of two, which has been encountered in similar context 
before [5TH55] . For A = 4 the temperature gives twice the Hawking temperature, i.e., 
T = 2 Th = 2 ( 87 tM) _ 1 . The literature is divided into two camps when it comes to the 
value of the temperature of radiation for a particle that tunnels through the horizon. 
Usually, those who favour T = 2 Th also favour the idea of dividing the entropy by 2 in 
order to satisfy Hawking’s original first law [53l;55|. However, according to Hawking’s 
original first law [50], for a static black hole [56] , 


Energy 

2 


Temperature x Entropy. 


( 66 ) 


Thus, if A = 4 then one should not divide the original entropy expression by 2 in order 
to get the correct internal energy on the l.li.s. of (j 66 ]h 
Also it is easy to check (159]) . 


s? = 



(67) 


This states that the system can be assigned a single temperature value which is a 
worldsheet-constant. One also finds the corresponding work (]46l) term and Gibbs free 
energy (ED term 


w = 



(47Tr 2 ) , 


Q = 


2(r — AM) 


167T 


( 47 rr 2 ) 


( 68 ) 


which are equal at the quasilocal thermodynamic equilibrium, i.e., 


Cj = = m. 


(69) 


4-2. Friedmann-Lemaitre-Roberts on- Walker (FLRW) geometry 
Now consider the FLRW metric in comoving coordinates 

a\t) 


ds~ = — dt + 


nr- 


dr 1 + a 2 (t)r 2 d£} 2 , 


(TO) 


where k = {— 1 , 0 , 1 } for open, flat and closed universes respectively and the comoving 
observer dyads are 


= ( 1 , 0 , 0 , 0 ), 


nU = 0 , 0,—,0 


ar 


E“ t = (O.AiZ'if 0,0 


n\= 0,0,0, 


ar sin 0 


(71) 


(72) 


Again note that the resultant thermodynamic potential densities do not change for 
observers with an arbitrary instantaneous radial boost with respect to the comoving 
observers. For such a set up, Helmholtz free energy density is 


f=\U a = m 2 


2 - 2 nr 2 - 2d 2 r 2 A 


a?r 2 


(73) 
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If we consider the equilibrium case where free energy takes its minimum value, one can 
find the equilibrium condition to be 


r VnTa 1 ' 


or 


r A = (ar) 


1 

s] H 2 + n/a 2,1 


(74) 


where H is the Hubble parameter. This corresponds to the location of the apparent 
horizon of the FLRW geometry [T811T9]. 

The internal energy density is found by isometrically embedding § into an Euclidean 
3-geometry and calculating its extrinsic curvature as k 0 = 2/(ar). Thus according to 
(1551) and (1551) . 


U = ( TS = 


2^H 2 + n/a 2 ' 


(75) 


with T = A/( 87 rar) = Aa J H 2 + n/a 2 / ( 87 t), and S = Area( 8 )/A = Air/ [A (H 2 + n/a 2 )]. 

For A = 4, this result matches the one of [151 [57H55] . where the temperature 
attributed to the apparent horizon is found to be T 4 = 1 / (27rr^). 

This temperature, assigned to the whole system, can be shown to be a worldsheet- 
constant by the variation 


^=2 \/ V 'VW 


A \ ~ ( A 
V“ 


87rar 


= 0. 


(76) 


When this condition holds, the work term (1461) and the Gibbs free energy (l5Tf are also 
found to be 


TF = 


2k + 2 ad + 2d 2 


16n 


( 47 t a 2 r 2 ) , Q = 


167T 


2d 

— + 


a r 2 a 2 


( 47 t a 2 r 2 ) . (77) 


By (ED, 


g = w = 


k + aa + a 2 


a z 


H 2 + n/a? 


(78) 


so that when the Friedmann equations are inserted we obtain 


g = 7F = 


1 - 3 p/p 


A.TT p 

3 


(79) 


where p is the pressure, p is the energy density of the perfect fluid and their ratio is 
u = p/p. Alternatively, we can compare the work required to create a quasilocal 2 - 
surface that encloses a system hlled with either vacuum energy (a; = — 1 ), stiff matter 
(u> — 1), dust (lo = 0) or radiation ( u> = 1/3). For the same value of the perfect fluid 
energy density, at equilibrium, the results state 


acuum ^ ‘l^Radiation 0 ) 


( 80 ) 
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meaning that a system filled with stiff matter has a greater tendency to store the 
potential relative work than a system filled with dust or radiation. Note that the 
surface tension is independent of the spatial size of the system in a FLRW spacetime, 
consistent with the fact that the FLRW geometry models a homogeneous universe. To 
see the differences with an inhomogeneous universe one may consider the Lemaitre- 
Tolman spacetime. 


4-3. Lemaitre-Tolman (LT) geometry 

The LT metric can be written in the comoving coordinates as 

R' 2 (t : r ) 


ds = —dt + 


1 + 2/1 (r) 


+ R (t, r)dQ , 


(81) 


where X' = and X — Qjf- for an arbitrary function X. One can choose a comoving 


^ and * - dx 
observer dual dyad 

ZW =(1,0,0,0), 




n \ = ( 0 , 0 ,-, 0 ), 


nO = 10,0,0, . 

d ' R sm 6 


(82) 

(83) 


Then the Helmholtz free energy density for such an observer is given by 


/ = 


1 


2 + 4 K - 2 R 2 
R 2 


(84) 


Minimizing R to obtain the condition for quasilocal thermodynamic equilibrium we find 


R 2 = 1 + 2/1, 


(85) 


which again gives the location of the generalized apparent horizon of the LT geometry 
|60j . In the absence of the cosmological constant, the evolution equation for the LT 
spacetime may be written as 

k2 = 2K+ If, ( 86 ) 

where M = M(r) is an arbitrary function which is said to play the role of the active 
gravitational mass within a constant radius shell in LT solutions. Therefore, another 
way of defining the apparent horizon is R(t,r ) = 2 M(r). Then after computing the 
internal energy density as k 0 = 2/R by (1391) and (|53]l . the internal energy becomes 

U = n:s = M{r), (87) 

with T = A /(87 tR) = A /(167 tM), and S = Area(E>)/\ = (16 ttM 2 ) /A. If we take A = 4, 
then the temperature assigned to the system takes the same value as the temperature 
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attributed to the apparent horizon in [BT]. The work term (HHl) and the Gibbs free 
energy (l5Tf are also found as: 


= 


1 

167T 


2 K' | 2 R | 2 RR> 
RR' + ~R + RRf 


(4vr R 2 ) , 


( 88 ) 


Q = 



2R 2 | 2(2/1+1) 

“R2" + ^ 


^2/1' - 2RR') 
R/L 


(4vr/? 2 ) . 


(89) 


Substituting (1861) and R(t, r) = 2 M(r) into the equations above gives 

g = %? = M(r)/V 2. (90) 

In contrast to the homogeneous cosmology, the surface tension in (188|) depends on the 
radial position of the quasilocal observers who define the inhomogeneous system. 

One can check the thermodynamic potentials of the LT system reduce to those of 
the FLRW and Schwarzschild spacetimes in the appropriate limit. In particular, for [60] 

R = a{t)r and M — J AnpR 2 R'dr (91) 

the relative work term (1901) agrees with the dust case of (]79l) for the FLRW geometry as 
expected. Likewise, if we take R = r and use the spatial derivative of the evolution 
equation (]86li to eliminate K'{r ) then at equilibrium (j8SJ) agrees with Q6Hj) for the 
Schwarzschild geometry. However, the general relative work (1901) of LT differs from 
the Schwarzschild case (169|) . on account of the competing terms in (1881) or (l89]h Recall 
that the Gibbs free energy is a measure of how much potential the system possesses to 
do work. In the static limit the second and third terms inside the square root in (l88|i 
vanish. In that case, the system stores all of the gravitational energy due to the spatial 
term —2 K'/ ( RR ') as potential work without being compelled to expend some of this 
energy as the system evolves in time. 


5. Local versus quasilocal equilibrium 

In Newtonian physics a system composed of particles which are in local thermodynamic 
equilibrium with each other is not always expected to be in global thermodynamic 
equilibrium. In the case of global equilibrium, one can assign single values of 
thermodynamic variables to the whole system. Likewise, in our case we have defined 
a quasilocal thermodynamic equilibrium condition so that one can assign a single 
temperature value to the whole system. 

In general, there is no reason for a system in local (or global) equilibrium to be 
in hydrodynamic equilibrium as well. One requires additional conditions for them to 
coincide [62], so we should be careful to distinguish these concepts. 
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A typical example of the local thermodynamics of matter fields on a curved 
background is given by Gao [63] who generalizes early work of Sorkin, Wald and 
Zhang ra to a generic perfect fluid. He investigates the connection between local 
thermodynamic equilibrium and hydrostatic equilibrium. Gao considers a collection of 
monatomic ideal gas particles with p = p(T), p = p(T) and s = s(p,n ) where s is 
the locally defined entropy density and n is the number density of the particles. By 
maximizing the local matter entropy, an equation for local hydrostatic equilibrium is 
obtained. This is not a unique way of defining the local thermodynamic equilibrium but 
this specific thermodynamic equilibrium condition also satisfies the local hydrostatic 
equilibrium. Also it is important to note that the fluid particles are not necessarily in 
a local thermal equilibrium here. 

Now let us consider the analogous problem of the conditions under which quasilocal 
thermodynamic equilibrium and quasilocal hydrodynamic equilibrium coincide for a 
general system containing both matter and gravitational energy contributions. We will 
assume that locally defined condition given by Green, Schiffrin and Wald [62] for matter 
fields also holds in the quasilocal case. This requires the quasilocally defined entropy to 
have its extremum value with respect to a total variation defined by (j58lh In order to 
compute this one can consider a generic spherically symmetric spacetime metric 

ds 2 = — A 2 (r , t)dt 2 + B 2 (r , t)dr 2 + i? 2 (r, t)<A2 2 , (92) 


where A(r, t) and B(r, t ) are arbitrary functions, R(r, t) is the areal radius of §, and our 
double dyad is 


^6 = 


A(r, t ) 


, 0 , 0,0 , 






(93) 

(94) 


Then, quasilocal Helmholtz free energy takes its minimum value when A 2 R' 2 = B 2 R 2 , 
and the total variation of the quasilocally defined entropy, fl54|) . at equilibrium is 


SS = -\/X7 a SV-S 


rf*E\d„ (TpJ= 0, 


= - a 


showing that the entropy is an extremum. This allows us to conclude that the quasilocal 
thermodynamic equilibrium and quasilocal hydrodynamic equilibrium should coincide. 
The interpretation of this result is crucial for the next section. 


6. Quasilocal virial relation 

Here we will sketch how the formalism above might be adapted to give a quasilocal virial 
condition which differs in character from previous attempts to define a virial theorem in 
general relativity [651168] . In previous studies only matter fields have been investigated, 
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in which the central object is the energy-momentum tensor defined at each spacetime 
point. However, a full description of the virial theorem in general relativity should also 
include gravitational energy, which cannot be defined at a point due to the equivalence 
principle. 

Ordinarily in classical mechanics the virial theorem is obtained by considering a 
system with motions confined to a finite region of space. If the potential energy is 
a homogeneous function of the coordinates, then the virial theorem gives a relation 
between the time averaged values of the total kinetic and potential energies ||69|. For 
such a system one can define a virial function G by [70] G(t) = '}Q l p l -r l , where fj are 
the coordinates and the p) are the momenta of the particles in the system. If G(t) is a 
bounded function, then the mean value of its time derivative is zero, i.e., 

) = ((m) ) + ^ ) = °- ( 95 ) 

Therefore, for a system under gravitational potential, the virial theorem reads 

2 (K.E) = - (. P.E) , (96) 

where (K.E) and (P.E) are the time averaged values of the total kinetic and potential 
energies. 

When it comes to including relativistic effects, however, this result changes. For 
the ultrarelativistic limit (v —>- c), the virial theorem takes the form 

{K.E) = - {P.E). (97) 

In astrophysics, the virial theorem is used when the thermal and the gravitational 
forces acting on an isolated system balance each other so that the system neither expands 
nor contracts. This state of the system is defined by its hydrostatic equilibrium [72]. In 
general the system is assumed to be composed of ideal gas particles which are in local 
thermal equilibrium with each other, which guarantees stability [62]. The value of the 
internal energy of the system, E in , is then equal to the ensemble average of the kinetic 
energies of the particles creating the system [73], i.e., 

E in = IGE, (98) 

where the overbar denotes the ensemble average at a given time. Note that the 
equipartition theorem is an application of the virial theorem. If thermal equilibrium 
coincides with hydrostatic equilibrium, then the temporal average of the kinetic energy 
of the total system becomes equal to the ensemble average of the kinetic energies of the 
particles at a given time (73], i.e., (K.E) = K.E |pp Consequently, for this case, 

one can rewrite (196]) and (1971) as 


2 E in = - 


(P.E) ( nonrelativistic ) 


(99) 
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E in — _ (P'E) ( ultrarelativistic ) (100) 

In onr case, quasilocal thermodynamic equilibrium is set by the minimum of the 
Helmholtz free energy which holds for systems with worldsheet-constant temperature 
and volume. This occurs when the mean extrinsic curvature of § is zero. Thus if one 
perturbs § along T, it neither expands nor contracts. Furthermore, the quasilocally 
defined entropy of the system then takes its extremum value. We can then expect the 
quasilocal thermodynamic and hydrodynamic equilibria to coincide. 

By analogy to the case of matter fields only, a virial theorem might hold also in our 
quasilocal formalism. In particular, we will suggest a natural bound, 

Ematter+grav EkLY E B ( 101 ) 

from which a quasilocal virial relation follows. Here E ma tter+grav represents the combined 
energy of the matter and gravitational fields, while E™ n is the minimum binding 
energy of the system. Furthermore, the interpretation of Ekly as an invariant proper 
mass-energy for generic spacetimes is agreed upon by many authors tHl[T2][36j. We 
demonstrate that (1 101 1) agrees with known results in particular limits, but will not 
attempt a formal proof of this bound or consequently of a virial theorem, which would 
require a detailed definition of E matter+grav in terms of suitable ensemble averages. 

We first recall that Bizon, Malec and 6 Murchadha m introduced a mass bound 
for a collection of spherical shells under collapse, given by 


M < M p - E b , (102) 

where M is the total energy of the shells, M p is the total proper mass and E B is the 
binding energy. The equality holds when binding energy is minimum, which turns out 
to be the Newtonian limit. Later, Yu and Caldwell [76] included this argument in their 
calculation of the binding energy of a Schwarzschild black hole and showed that 


M — M p — E% in , 


(103) 


and M p = E B y in the Schwarzschild geometry for static observers at any r. Since 
E B y = Ekly for static observers in the Schwarzschild geometry for observers who 
are instantaneously at rest, (1103|1 is seen to coincide with (HOip in the Schwarzschild 
geometry. 

Now let us specialize to observers at quasilocal hydrostatic equilibrium in the 
Schwarzschild geometry, at r = 2M, where the Kijowski-Liu-Yau energy gives 


Ekly — 



2 \ 1 — 


2 M 
r 


r 



= 2 M. 


(104) 
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Since the internal energy, E’p (J7J), at equilibrium corresponds to the usable matter + 
gravitational energy, the l.li.s. of (1101(1 should read 


zp _ jpnn 

^ matter -\-grav ^ K 




= M. 


Hence Eff m = M. The minimum binding energy of such a system can also be found by 
calculating the work done in bringing the spherical mass shells from infinity to r = 2 M 
by observers who are instantaneously at rest in the Newtonian limit ra. Thus at 
hydrostatic equilibrium, 


Tp j_ Tpmin 

^ matter-\-grav ^ B ? 


(105) 


which is analogous to the virial relation (HOOD in the ultrarelativistic limit, since binding 
energy is negative of the potential energy, and only the so-called reference term survives 
in E i p at quasilocal thermodynamic equilibrium. 

Now we will generalize this result by assuming that (1101(1 holds for generic 
spherically symmetric spacetimes at quasilocal hydrodynamic equilibrium. At quasilocal 
hydrodynamic equilibrium yj2 (J^) = pk 2 — l 2 = 0 , by (JTJ) and (J8]), 


E, 


matter-\-grav 


rpin j_ 

- e k - 


167T 


Hence from IIOil we find 



77 ’ 77 ' Tpmin J_ 

t- J KLY -E- J matter-\-grav *- J B 


167T 


(106) 

(107) 


(108) 


where k 0 = 2/R(r, t ) is the extrinsic curvature of § when embedded in Euclidean 3-space 
and R(r, t ) is the areal radius of S. Therefore, by (1106(1 and (11081) one obtains (11051) as 
a virial relation for any spherically symmetric distribution at quasilocal hydrodynamic 
equilibrium whose matter and gravitational contributions to the total content cannot 
be decoupled. 

A key inference is that the proper mass-energy, usable matter-energy and the 
binding energy of a system make most sense when referred to measurements made by 
the same set of quasilocal observers. Therefore the idea of comparing certain energy 
definitions for systems with different sizes in a given spacetime can lead to paradoxes. 
For example, Frauendiener and Szabados argued that .if the quasi-local mass 

(Ej kly) should really tend to the ADM mass as a strictly decreasing set function near 
spatial infinity, then the Schwarzschild example shows that the quasi-local mass at the 
event horizon cannot be expected to be the irreducible mass. ” This is true simply because 
a system with a spatial size coinciding with the event horizon has different binding 
energy requirements than the one whose spatial size tends to infinity. The latter one 
























Quasilocal energy and thermodynamic equilibrium conditions 


23 


has zero binding energy, because no work has to be done by the system defined by 
quasilocal observers located already at infinity. The authors continue with the statement 
“...there would have to be a closed 2-surface between the horizon and the spatial infinity 
on which the quasi-local mass would take its maximal value. However, it does not seem 
why such a (geometrically, and hence, physically) distinguished 2-surface should exist.” 
Here we note that such a closed spatial 2-surface does exist with a location matching 
the apparent horizon. It encloses a system whose quasilocal thermodynamic equilibrium 
coincides with quasilocal hydrodynamic equilibrium. This might also serve as a physical 
interpretation for a generalized apparent horizon for the case when its location matches 
the one of a marginally outer trapped surface. The outgoing null rays of a system 
enclosed by such a trapped surface do not tend to leave the system because systems in 
quasilocal thermal equilibrium simply do not radiate. 

7. Discussion 

When the equilibrium thermodynamics of horizons was first introduced in the 1970s 
M, the quasilocal energy definitions that we have today were unknown. It is now 
known that the physically relevant boundary Hamiltonian of general relativity lies on a 
closed 2-dimensional spacelike surface, S, of a spacetime domain [101114] , which we call 
the screen. In this paper we have focused on a spherically symmetric system enclosed by 
a screen, S, as the central object of gravitational thermodynamics rather than horizons. 

Isolated systems are natural objects in classical thermodynamics. In general 
relativity, however, no system can be totally decoupled from the rest of the universe due 
to the nonlinear nature of the gravitational interaction. The systems we consider in this 
paper have arbitrary size and are generally in nonequilibrium with their surroundings. 
Only after quasilocal thermodynamic equilibrium conditions are introduced does it 
follow that the screen is located at the apparent horizon of [15], where the standard 
equilibrium thermodynamic laws apply. 

We believe that this approach may ultimately prove useful in general relativity, 
since the issues associated with quasilocal gravitational energy on one hand, or with 
gravitational entropy on the other, are generally studied in isolation. In fact, the problem 
of gravitational entropy is so complex that often researchers simply seek definitions in 
terms of geometric quantities which are nondecreasing with time [78b8(Tj . giving rise to 
a “second law”, without directly investigating whether the entropies so defined obey 
any of the other properties one might reasonably demand of a genuine thermodynamic 
potential. The fact that the second law of classical thermodynamics can be viewed as 
a consequence of entropy not being rigorously defined in nonequilibrium [40] is usually 
overlooked. 

On account of the equivalence principle, statistical macroscopic properties of the 
gravitational field are necessarily nonlocal. To interpret quasilocal gravitational energy 
in terms of thermodynamic laws it is necessary to have a measure of the “work done” by 
the tidal forces on the screen associated with the quasilocal observers. For this reason, 
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we have adapted the generalized Raychaudhuri equation of an arbitrary dimensional 
worldsheet embedded in an arbitrary dimensional spacetime [13], to the special case of 
a 2-dimensional timelike surface, T, (orthogonal to § at every point), which we embed 
directly into 4-dimensional spacetime. 

The mean extrinsic curvature of S, that appears in the quasilocal energy definitions 
[T0TFT2] . gives the expansion of 8 when it is perturbed along T. Degrees of freedom 
living on T are , therefore, understood to be those which describe the changes of 
thermodynamic potentials, while the degrees of freedom living on the screen, 8, are 
required to consistently define a system. Hence, in order to write the first law, the total 
variations of the thermodynamic variables are taken along the 2-surface T rather than 
variations along the integral curve of a single vector. 

It is known that quasilocal energy definitions which involve the mean extrinsic 
curvature of 8 are invariant under radial boosts. In our formalism this boost invariance 
holds for all thermodynamic potentials that appear in the generalized Raychaudhuri 
equation (l27lh once the definitions (l33lh (1451) and (l50li are made. This is possible on 
account of spherical symmetry. 

Spacetime geometry is the source of gravitational energy while also defining the 
relationships between the clocks and rulers that must be used to investigate the 
system. One might wonder how quasilocal thermodynamic relations could be treated in 
spacetimes more general than spherically symmetric ones. In the absence of symmetries 
there are no unique definitions of energy in general relativity. Conservation laws, by 
definition, require that some property does not change. Therefore only with some 
symmetry condition can one define a system consistently by keeping certain degrees 
of freedom constant. 

Since gravitational energy and related macroscopic variables are not local, whether 
any particular definition has utility depends on the consistent definition of 2-surfaces, 
8, which are best adapted to the spacetime in question. Our spherically symmetric 
formalism might be easily extended to situations which are approximately spherically 
symmetric, in a perturbative scheme. Furthermore, we believe that similar reasoning 
to ours could also apply to spacetimes with other symmetries, such as axial symmetry. 
In that case one should be able to introduce additional quantities, which account for 
rotational energy, for example. I 11 such a case the first and third terms on the r.h.s. 
of the generalized Raychaudhuri equation (1271) are nonzero, making a thermodynamic 
interpretation considerably more complicated, and is left to future investigations. 

As an application of our formalism, we have sketched a natural bound involving the 
quasilocal gravitational energy plus matter fields, which might suggest a virial relation. 
To rigorously prove a virial theorem requires that we have a proper understanding of 
the degenerate states of matter and gravitational fields contained within the screen, 8, 
which are consistent with the same worldsheet-constants on 8. Such an understanding 
of course requires going far beyond the present paper, as it effectively means probing 
fundamental questions related to the holographic interpretation which are important to 
both statistical and quantum gravity. 
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Other questions for future work relate to the question of nonequilibrium quasilocal 
gravitational thermodynamics. The Helmholz free energy, (I35[) . is defined for all states 
of the system, whereas the other thermodynamic potentials have only been defined at 
quasilocal thermodynamic equilibrium. This is simply because equilibrium was defined 
by minimization of jF, which may not be the only way to define a useful quasilocal 
thermodynamic equilibrium condition. Other types of equilibrium conditions could be 
applied to the generalized Raychaudhuri equation. 

For the case of thermodynamic nonequilibrium, the Raychaudhuri equation of the 
worldsheet (1271) should still quantify the energy fluxes into and out of the system. 
However, in that case, the existence of thermodynamic variables is not guaranteed and 
their consistent definition becomes murky even in classical thermodynamics. Losing 
the linearity condition among the thermodynamic variables makes their interpretations 
much more difficult and one might need additional geometrical structures to investigate 
processes involving their generation and dissipation. 

Recently, Freidel [8] presented an approach to study nonequilibrium thermodynam¬ 
ics by using geometrical objects in a 2+2 formalism. Our language is similar to his in 
terms of the quasilocal nature of the thermodynamic potential densities introduced on 
a screen S. However, our formalism differs fundamentally in character by the existence 
of a worldsheet T on which both of the degrees of freedom are treated equally in terms 
of their roles in evolving the potentials. Whether or not investigation of this difference 
provides a passage from quasilocal equilibrium to nonequilibrium thermodynamics is a 
point of interest. 

In conclusion, we have introduced a new formalism to link two areas which are 
often treated separately - quasilocal energy and horizon thermodynamics - by adapting 
a geometrical approach of Capovilla and Guven [T3] to define a quasilocal equilibrium 
condition. We have investigated specific examples and potential applications in the case 
of spherical symmetry. This already sheds fresh light on some paradoxes relating to 
quasilocal energy. Although our paper is by nature exploratory, we believe that the 
steps taken here can be further built upon to potentially tackle important questions 
that generally are not currently considered, but should be, in gravitational physics. 
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